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ABSTRACT 

The development of a unifying framework among direct 
numerical simulations, large-eddy simulations, and statisti- 
cally averaged formulations of the Navier-Stokes equations, 
is of current interest. Toward that goal, the properties of 
the residual (subgrid-scale) stress of the temporally filtered 
Navier-Stokes equations are carefully examined. Causal 
time-domain filters, parameterized by a temporal filter width 
0 < A < oo, are considered. For several reasons, the 
differential forms of such filters are preferred to their corre- 
sponding integral forms; among these, storage requirements 
for differential forms are typically much less than for inte- 
gral forms and, for some filters, are independent of A. The 
behavior of the residual stress in the limits of both vanishing 
and infinite filter widths is examined. It is shown analyti- 
cally that, in the limit A — y 0, the residual stress vanishes, in 
which case the Navier-Stokes equations are recovered from 
the temporally filtered equations. Alternately, in the limit 
A — y oo, the residual stress is equivalent to the long-time 
averaged stress, and the Reynolds-averaged Navier-Stokes 
equations are recovered from the temporally filtered equa- 
tions. The predicted behavior at the asymptotic limits of 
filter width is further validated by numerical simulations of 
the temporally filtered forced, viscous Burger’s equation. Fi- 
nally, finite filter widths are also considered, and a priori 
analyses of temporal similarity and temporal approximate 
deconvolution models of the residual stress are conducted. 

INTRODUCTION 


The Navier-Stokes equations can be solved numerically 
to predict turbulent flows; however, due to the enormous 
computational expense required to extract a solution from 
these equations for flows of engineering interest, it has 
been necessary in most cases to revert to alternate formula- 
tions. For current purposes, three computational approaches 
are considered: direct numerical simulation (DNS), large- 
eddy simulation (LES), and Reynolds-averaged Navier- 
Stokes (RANS) computations. These differ primarily in the 
level of approximation required to achieve closure. 

For LES, the separation of the field variables into re- 
solved and unresolved (spatial) scales is effected by filtering 
the fields with a low-pass filter. Filtering the momentum 
equations generates residual (subgrid-scale) stresses that re- 
quire closure either by modeling or approximation. Recent 
advances such as dynamic modeling (Germano, 1991) and 
deconvolution methods (Domaradzki and Saiki 1997, Stolz 
and Adams 1999) have made LES practical for application 
to certain flows of engineering interest (Moin and Jimenez, 
1993). 

Long-time averaging of the Navier-Stokes equations re- 
sults in the RANS equations for the time-independent mean 
state. RANS methodology is generally applied to statisti- 
cally steady flows. To close the system of equations, a model 
is needed for the Reynolds-stress tensor. Although RANS is 
computationally appealing, it places a heavy burden on the 
Reynolds-stress model, which must incorporate the effects 
of all the unsteady motions upon the mean. 

While the formal linkage of the LES and RANS equations 



has been well established (Germano 1999, 2001), it is of in- 
terest to investigate whether this linkage can be extended 
practically by developing filtering and averaging procedures 
that yield mutually consistent solution fields. A possible 
unifying context for these methodologies is afforded by fil- 
ter theory. However, the linkage between LES and RANS 
may be more natural within the context of time-domain fil- 
tering rather than the traditional spatial filtering commonly 
used in LES. Accordingly, the present study focuses on the 
temporally filtered Navier-Stokes (TFNS) equations and the 
resultant residual-stress fields. 

In the next section causal time-domain filters are dis- 
cussed, and differential forms are derived for two candidate 
filters: an exponential filter and a Heaviside filter. The 
TFNS equations are formulated and analyses of the asymp- 
totic behaviors of the residual stress for limiting values of 
filter width are presented. Finally, temporal residual-stress 
models are proposed for the case of finite filter width. Fi- 
nally, the numerical solution of the forced, viscous Burger’s 
equation is used to further validate the analytical results 
as well as to evaluate the proposed temporal residual-stress 
models by a priori analyses. 


TEMPORALLY FILTERED NAVIER-STOKES EQUATIONS 

Time-domain filters are classified as causal or acausal 
depending upon whether they are applicable to real-time or 
a posteriori data processing, respectively. The interest here 
lies in real-time applications for which only causal filtering is 
appropriate; accordingly, the focus in this study is restricted 
to causal filters. While aspects of time-domain filters have 
been discussed previously in this context (Pruett 2000), it is 
worthwhile to reiterate some fundamental relationships for 
completeness. 

Let f(t) be a continuous function of time t. A causal 
linear filter is readily constructed by the integral operator 

/( t ; A )= f G(t - t;A)f(r)dT, (1) 

J — CO 

where G is a parameterized filter kernel, and the parameter 
A is the filter width. In general, admissible kernels must 
satisfy the following property: 


G ( t; A ) = A/ 


t 

A 


(2) 


where g is any integrable function such that 


g(t) > 0, 


g(t)dt = 1 


— CO 


and p(0) = 1. 


( 3 ) 


The non-negativity and normalization constraints in Eq. (3) 
imply that 

lim g(t) = 0, (4) 

t—¥ — o o 

and suffice for G to approach a Dirac delta function as its 
parameter A — y 0; that is, 


lim f(t; A) 
A-HD 
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A-HD 
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t;A)f(r)dT 


f S(r — t)f (r )dr 
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( 5 ) 


Two examples of simple, useful filters are obtained by use of 
an exponential function and a Heaviside function as kernels. 


For the exponential function, the kernel is 

g(t) = exp (t) -»■ G(t; A) = ex P 0/ A ) , 
and the resulting integral operator in Eq. (1) is 

/( t;A ) = X f eXp (~A~) ^ T ' )dT ' 

J — CO 

If g(t) = h{t + 1), with h the Heaviside function, then 

/( t;A ) = 3 - f f(r)dT. 

^ Jt - A 


(6) 

( 7 ) 

(8) 


A drawback of the integral formulations just presented 
is the need to retain the long-time history of the solution 
field. However, by considering instead differential forms of 
the filter operators, storage requirements are reduced signifi- 
cantly, subject to the intrinsic storage needs of the numerical 
time-advancement scheme itself (for example, low-storage 
Runge-Kutta). By differentiating Eqs. (7) and (8), the dif- 
ferential forms of the exponential and Heaviside filters are 
given by 


/(*; A ) = 
/(*; A ) 


/(*) - /(*; A ) 
A 


and 


A 


( 9 ) 

( 10 ) 


respectively. 

When causal filtering is applied to a temporally dis- 
cretized problem with a time increment of At, the action 
of the filter is naturally parameterized by the filter-width 
ratio r defined as 


In order to illustrate the discrete differential filtering pro- 
cess, a 27r-periodic time series is processed by the exponential 
differential filter given in Eq. (9). The time series is gen- 
erated from a -3/2 power-law decay in Fourier frequency 
space, and the phases are assigned randomly. The continu- 
ous signal is then sampled at a rate of 512 per period and 
replicated for three periods. The filtered time series, /, is 
then generated by solving Eq. (9) from the initial condi- 
tion /( 0; A) = /( 0). There are many appropriate numerical 
integration schemes. Because the right-hand side of the dif- 
ferential form of a linear filter is itself linear, fully implicit 
Adams-Moulton methods are particularly attractive because 
of their accuracy and efficiency. Here, standard fourth- 
order Adams-Moulton methodology is used. The method 
is started with initial steps of orders one, two, and three, 
respectively. 

The filter-width ratio, r, is the only parameter of the 
differential filter. In general, the larger the value of r, the 
more dissipative the filter. (In this context, a “dissipative” 
low-pass filter is one with significant and broad-band atten- 
uation of high-frequency Fourier harmonics.) The method 
remains viable for all values of filter-width ratio (0 < r). 
For r 0, the evolution equation becomes stiff, and small 
time steps are necessary. Figure 1 compares the filtered time 
series with the unfiltered signal for selected values of the 
filter-width ratio r. As r increases, the output time series 
becomes smoother and its amplitude diminishes due to the 
removal of energy at the higher frequencies. As is typical 
for causal filters, high levels of numerical dissipation gener- 
ate significant phase lag in the output relative to the input. 
Figure 2 compares the input signal with the original out- 
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be independent of the choice of the filter. To make explicit 
its formal dependence upon the specific temporal filter, the 
residual stress is denoted by [tr]^(A) where appropriate. 


Limiting Behaviors 

Of interest is the effect of filter width A on the residual 
stress [tr]*j. It is easily shown that the [tr] vanishes in 
the limit A — y 0. In this limit, the kernel function reduces 
to a Dirac delta function (see Eq. (5)) so that 

lim [trL ? ' (t, x'. A) = lim (ui u n — UiUj) 

A— H3 L A-H) J 

= lim (uiUj — uiu-j) = 0. (15) 
A-H) 


Figure 1: Differentially filtered time series /(i). 

put signal and with an output that is phase compensated 
by r time steps. The phase-compensated signal is an ex- 
cellent representation of the input, minus its high-frequency 
components. 



The vanishing of the temporal residual stress, coupled with 
the replacement of the other filtered quantities by their unfil- 
tered counterparts, leads to the recovery of the Navier-Stokes 
equations from the TFNS equations in the limit A — y 0. 

The other limit of interest is A — y oo. However, before 
examining the behavior of the residual stress in this limit, 
it is useful to examine some characteristics of the filtered 
velocity field itself. It follows from the differential forms of 
either the exponential or Heaviside differential filters given 
in Eqs. (9) or (10) that 

du,i 

lim -(i, x; A) = 0, (16) 

A — >■ co dt 

where both Ui and u{ are assumed bounded. The condition 
above establishes that ui(t,x; oo) is actually independent of 
time t. (In fact, for Eq. (16) to hold, it suffices that Ui 
is bounded and that | < g |/ (t)| is integrable on ( — oo,0]. Thus, 
Eq. (16) applies to a wide class of filters.) For the causal 
temporal filter defined in Eq. (l) with the Heaviside kernel 
(for convenience), v,i(t,x; oo) can be written as 


Figure 2: Original, exponentially filtered, and phase-com- 
pensated time series f(t) for r = 32. 

Temporal, causal filtering of the Navier-Stokes equations 
using Eq. (l) leads to the following form of the TFNS equa- 
tions: 


:du L 

dx 3 



d(u t uj) 

dxj 


o (12) 

_d L + v d 2 ~ t _£M^ (13) 

dxi dxjdxj dxj 


where Ui is the velocity, p is the pressure, and v is the kine- 
matic viscosity. An overbar denotes a temporal grid-filtered 
quantity, and [tr] represents the temporal residual-stress 
tensor defined as 


lim Ui(t,x; A) 
A— ^oo 


i (0, x', oo) = lim 
A— fr-cx 


1 

A 


i (t, x)dr 


(17) 

Equation (17) holds for any filter for which 77(0) = 1, which 
is typical of low-pass filters. The right-hand-side of Eq. (17) 
simply defines the long-time average of the variable Ui (t, a;), 
which, for a stationary process, is equivalent to the ensemble 
average according to the ergodic hypothesis. That is, for a 
stationary flow 


Ui (0, x; oo) = E\ui (i, a;)}, (18) 

where 77{ } denotes the expected value (or ensemble av- 
erage). However, Eq. (16) has shown that v,i(t, x;oo) is 
constant with respect to time so that 


[trJjj = it; it j - UiUj. ( 14 ) 

Provided that filtering and differentiation operations com- 
mute, the TFNS equations are formally identical to the 
spatially filtered Navier-Stokes equations. As pointed out 
previously by Pruett (2000), commutativity is natural for 
temporal filters but remains problematic for spatial ones 
(Blaisdell 1997, Vasilyev et al. 1998). It is now recognized 
(Pruett 2001, Pruett et al. 2001, De Stefano and Vasilyev 
2001) that this formal equivalence does not imply quantita- 
tive equivalence of the residual-stress fields. In general, for 
spatial or temporal grid filters, the residual stress depends 
strongly upon the filter, particularly upon its filter width and 
order property, which influence both the magnitude and the 
distribution of the residual stress. The implication of this 
growing awareness is that the residual-stress model cannot 


Ui{t,x; oo) = i4*(0,a;; oo) = E{ui(t,x)f = Ui{x), (19) 

and 

lim — ui (t, x; oo) = — E{ui(t,x)} = 0. (20) 

A— »-oo dt dt 

In the current time- filtered approach, Eq. (19) provides the 
link between the resolved motions of the variable Ui(t,x) 
and the ensemble mean U,( x). Because the variable Ui(t,x) 
can be partitioned either into a sum of resolved Ui(t,x', oo) 
and temporally unresolved motions Ui(t,x), or into a sum 
of time mean U,( x) and fluctuating u^{t,x) quantities, it 
follows from Eq. (19) that 

u,i (t, x; oo) = u\ (t, x). (21) 

In addition to the equality between the resolved and mean 
fields in the limit, Eq. (21) shows the linkage between the 




temporally unresolved and fluctuating motions. With these 
results, it is now possible to examine the limiting behavior 
of the residual stress. 

By the linearity of the filter operator, the residual stress 
defined in Eq. (14) can be written as 

lim [r rJjj (A) = limA^oo \(u t Uj + u t Uj 
A — )■ co L 

-\-u t Uj + u t uj) — u t Uj J , (22) 

where the instantaneous velocity field has been partitioned 
into resolved and unresolved parts. Because Eqs. (18) and 
(19) establish an equality between the resolved and ensemble 
mean fields, and the residual and fluctuating fields, respec- 
tively, Eq. (22) can be simplified to 

lim [t r ]^(A) = E{u^E{u{} + u' t E{u 3 } + u^uA 

A— too J J 

= E^uj} = rij. ( 23 ) 

That is, for a stationary flow the residual stress ([tr]^) 
asymptotically approaches the Reynolds stress (r t j) as A — y 


Finite Filter Width 

For finite filter width, the residual stress represents the 
dynamics of a broad spectral range of motions. The TFNS 
equations provide a governing set of equations suitable for 
time-filtered LES (Pruett, 2000) (or TLES), for which (ac- 
curate) modeling of [t r ]^ is required for closure. 

The temporal variants of two well-known residual-stress 
models for [tr]^ are considered: Bardina’s scale-similarity 
model (Bardina et al., 1980) (SSM);andthe approximate de- 
convolution model (ADM) of Stolz and Adams (1999). The 
time-filtered counterparts of these models are referred to as 
the temporal scale-similarity model (TSSM) and the tempo- 
ral approximate deconvolution model (TADM), respectively. 

Consider first a TSSM that is formally equivalent to the 
Bardina model (Bardina et al., 1980), 

~ u t Uj — u t Uj, (TSSM). (24) 

As in the Bardina model, the same (temporal) filter width 
is used for the primary and secondary (test) filters. Next, 
the TADM considered is formally equivalent to the second 
of the ADM models presented by Stolz and Adams (1999), 

[tr]*j % viVj - v t v 3 , (TADM), (25) 


the full TFNS equations, the wide range of parameter val- 
ues considered here renders such analyses cost prohibitive. 
However, it is possible to illustrate the dependence of the 
residual stress upon the temporal filter width, in general, 
and the asymptotic behaviors discussed previously, in par- 
ticular, by simulations of a spatially one-dimensional model 
problem. To this end, consider the forced, viscous Burger’s 
equation (VBE), written in the form 


du 1 d(u 2 ) 

dt 2 dx 


'FT + 

dx z 


(0 < x < 27t). (26) 


with u(t,x) a velocity, f(t,x) an imposed forcing function, 
and v a viscosity. The initial condition is -u(0, x) = 1 . With- 
out forcing, the initial condition results in a velocity field 
that is constant for all time and space. Moreover, any per- 
turbations of that field decay toward zero, so that constant 
“stirring” is required to maintain high-intensity fluctuations 
(Eswaran and Pope, 1988). This equation can be solved 
accurately by a Galerkin Fourier spectral method in space 
coupled with classical 4th-order Runge-Kutta time advance- 
ment. A Fourier ansatz is assumed for u and substituted 
into the governing equation. This results in a system of cou- 
pled ordinary differential equations for the complex Fourier 
coefficients £/*., k = — n/2, ..., — 1, 0, +1..., +n/2. (Due to 
conjugate symmetry, only nj 2 + 1 positive modes are solved 
for explicitly.) The equations are coupled through their non- 
linear terms, which are evaluated exactly in Fourier space by 
Cauchy products. Hence, explicit de-aliasing is unnecessary. 

For this forced case, n = 256, and each Fourier mode in 
the band 1 < k < kf is independently subjected to periodic 
forcing (t) such that Fk(t) = A4 >k exp with real fre- 
quency LUk = ku>. The band limit kf = 32, the fundamental 
frequency to = 1, and the amplitude A = 0.4 (the same for 
all modes) are input parameters, and the time increment 
is 0.005 throughout. The complex phases (f>k = exp(to'fc) 
are assigned initially by random numbers a ^ uniformly dis- 
tributed on [0, 2 tt] . Thereafter, they remain fixed. As will be 
shown, after a long-time evolution, a statistically steady flow 
results. Because, at small is, the viscous Burger’s equation 
admits solutions with steep shock fronts, only a moderately 
large value of is is practical. For the value is = 1/300 and 
forcing distribution, the flow is highly resolved in both time 
and space, with Fourier amplitudes at the highest wavenum- 
bers on the order of 10 — 10 . 

Causally filtering the forced VBE results in the follow- 
ing equation, which can be considered as a one-dimensional 
analog of the TFNS equation given in Eq. (13): 


where v t is an approximate deconvolution of u t ; that is, v t 
approximates u t based upon approximately defiltering (de- 
convolving) u t . Following Stolz and Adams, the zeroth- 
and first-order deconvolutions of u t yield v t = u t and 
Vi = 2 u t — u t , respectively. Higher-order (and more ac- 
curate) deconvolutions are possible. Note that the TADM 
(ADM) generalizes the TSSM (SSM), because the zeroth- 
order deconvolution is the TSSM. (Appropriately, Stolz and 
Adams (1999) refer to the second of their ADM models as 
the generalized SSM model.) 

In the next section, an a priori analysis of the predic- 
tive capability of these two residual-stress models will be 
performed. 


FORCED VISCOUS BURGER'S EQUATION 

While it is desirable and ultimately necessary to validate 
the analytical results previously established in simulations of 


du 1 d(uu) d 2 u - 1 9[tr] 

— T = is — ~ T f [t. x) — 

dt 2 dx dx 2 2 dx 


(0 < x < 27r). 

(27) 


Filtering results in the appearance of a residual stress given 

by 


[tr] = UU ~ UU. 


(28) 


In Fig. 3, the instantaneous unfiltered velocity field, ob- 
tained from the solution of Eq. (26) at t = 10 (At = 0.005), 
is compared with the filtered field, which satisfies Eq. (27), 
for selected values of the filter-width ratio r. Clearly, filter- 
ing in time to remove high frequencies effects the removal of 
energy at high wavenumbers as well. 

As implied previously, the behaviors of the residual stress 
for limiting values of the temporal filter width A are key 
results of the temporally filtered methodology being studied. 
To illustrate these predicted behaviors in the limits A — y 0 
and A — y oo, the model problem is particularly useful. 
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Figure 3: Instantaneous and causally filtered velocity fields 
at t = 10 for filter- width ratios r = l,r = 4, r = 8, r = 16, 
and r — 32. 
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Figure 4: Instantaneous residual stress [tr] at t = 10 for 
filter- width ratios r — 1, r = 8, and r — 32. 

Limiting Behavior of Residual Stress 

The behavior in the limit A — y 0 can be verified nu- 
merically by using successively smaller temporal grid-filter 
widths to process the numerical solution u(t,x) of the VBE. 
The exact residual stress (28) is evaluated to the accuracy 
of the numerical scheme by solving, in addition to Eq. (26), 
the filter evolution equations (cf. Eq. (9)) 


du 

dt 



and 


duu 

dt 


uu — uu 

A 


(29) 


from initial conditions ti(0,a:) = it(0,a;) and uu( 0, x ) = 
u 2 (0,x). Here, these equations (29) are advanced in time 
using the standard fourth-order Adams-Moulton method. 
(The fourth-order Runge-Kutta methodology used to ad- 
vance the VBE would also be suitable for all the filter 
equations; however, following the Runge-Kutta update of 
the solution by the 4th-order Adams-Moulton updates of the 
filtered quantities has the algorithmic advantage of compart- 
mentalizing the code.) While these Eqs. (29) apply to the 
exponential filter, an analogous set could be derived for the 
Heaviside filter. 

Figure 4 compares the exact, instantaneous residual 
stress [tr] at t = 10 determined from Eqs. (26) and (29) 
for selected values of the filter- width ratio r (= A/ At). As 
expected, the amplitude envelope of the residual stress tends 
toward zero as r decreases. 

The behavior of the velocity and residual stress fields in 
the limit of A — y oo can also be analyzed. As was shown 
in previously, the limiting form of the residual stress [tr] 
approaches the long-time average stress field t. While this 
can be shown by considering successively larger values of the 


Figure 5: Instantaneous residual stress [tr] at t = 240 for 
selected values of filter- width ratio r: long time stress 


filter- width ratio r, it is first necessary to establish that the 
solution of the forced, viscous Burger’s equation evolves to a 
statistically steady (stationary) state. Thus, it is necessary 
to verify the stationarity of the numerical solution, because 
the equality of the residual stress and Reynolds stress in the 
long-time limit is based on this assumption. (See Eqs. (19) 
and (20).) 

The long-time average and the spatial average of the in- 
stantaneous velocity field u(t, x ) are both equal to unity. The 
fluctuating field u'(t,x) is extracted at each time step sim- 
ply by subtracting this mean value from u(t, x ). It is verified 
that the (long-time) solution of the forced viscous Burger’s 
equation is indeed stationary. With the stationarity of the 
solution of the forced VBE established, it is possible to eval- 
uate the effect of large filter width on the residual stress. In 
Fig. 5 the long-time averaged stress r is compared with the 
instantaneous residual stress [tr](A) at t = 240 for different 
values of the filter- width ratio r. The stress r is computed 
by averaging over an interval of duration A = 20 during the 
period in which the flow is essentially stationary (t = 240 
to 260). As expected, the residual stress [tr], computed in 
real time using the exponential filter, appears to converge 
toward the value of r as A becomes large. 


Residual Stress for Finite Filter Width 

With the limiting behavior of the filter-width A on [tr] 
established for the forced VBE, it remains only to evalu- 
ate the behavior of [tr] for finite filter width. As described 
previously, the exact residual stress is extracted from the 
solutions of Eqs. (26) and (29). The modeled residual stress 
can be obtained from these equations by further appending 
the evolution equations 


du u — u 

duu uu — uu 

(30) 

dt 

A ’ 

dt 

A 

dv v — v 

dvv vv — vv 

(31) 

dt 

A ’ 

dt 

A 


Equation (30) is used in conjunction with the TSSM, sub- 
ject to the initial conditions ti(0, x ) = it(0, x ) and uu{ 0, x) = 
it 2 ( 0,a;). For the TADM, both sets, Eqs. (30) and Eq. 
(31) are involved, subject to the additional initial conditions 
•i;(0,a:) = it(0, a?) and tn7(0, x') = it 2 ( 0,a;). As before, these 
differential filter equations are advanced by the Adams- 
Moulton method. 

In Fig. 6, the exact ([tr]) and modeled ([tr](TSSM) and 
[tr](TADM)) residual stresses are compared at t = 20. At 
this instant, the flow statistics are still evolving in time, 







Figure 6: Instantaneous exact and modeled (TSSM, TADM) 
residual stresses ([tr]) at t = 20. 

for in the finite-A case, there is no reason to presuppose 
stationarity. A fairly dissipative filter of ratio r = 16 is 
used for the a priori analysis. Note that for both the exact 
and modeled residual stresses, [tr] > 0 at all times. This 
realizability property (Vreman et al., 1994) is a consequence 
of the positivity of the filter kernel established in Eq. (3). 

Both the TSSM and TADM correlate relatively well with 
[tr], with correlation coefficients on the order of 0.8 and 0.9, 
respectively. Correlations, however, reflect distribution but 
not amplitude. In general, the TADM has a higher corre- 
lation, and its amplitude tends to be more nearly correct. 
That the TADM performs well with only first-order decon- 
volution is surprising (as Stolz and Adams (1999) employ 
5th order). 

CONCLUSIONS 

The behavior of the residual stress of the temporally 
filtered Navier-Stokes (TFNS) equations was studied for a 
class of differential, causal time-domain filters parameterized 
by the temporal filter width A. The effect of filter width on 
the residual stress was examined for the asymptotic limits 
A — y 0 and A — y oo and for the case of finite filter width. 
It was shown analytically that, in the limit A — y 0, the 
residual stress vanishes so that the Navier-Stokes equations 
are recovered from the temporally filtered equations. Alter- 
nately, in the limit A — y oo, for a statistically steady flow, 
the residual stress asymptotically approaches the Reynolds 
stress, and the Reynolds-averaged Navier-Stokes equations 
are recovered from the temporally filtered equations. These 
asymptotic results were verified numerically through sim- 
ulations of the temporally filtered forced, viscous Burger’s 
equation. For the case of finite filter widths, two residual- 
stress models were considered that are temporal analogs of 
spatial SGS-stress models. These were a temporal scale sim- 
ilarity model (TSSM) and a temporal approximate deconvo- 
lution model (TADM). A priori analyses of these models 
were performed using highly accurate numerical solutions 
of the filtered forced, viscous Burger’s equation. The mod- 
els were found to approximately replicate the exact residual 
stress. 

It has been shown analytically that the residual stress 
of the TFNS equations is strongly dependent upon the tem- 
poral filter width. This fact, coupled with computational 
results from simulating the forced, viscous Burger’s equation 
over a wide range of temporal filter widths, suggests that full 
simulations of the TFNS equations should behave like DNS 
for small temporal filter widths and like RANS for very large 
ones. For finite filter widths the formulation describes a tem- 


porally filtered LES or TLES. These results have provided 
a bridging mechanism between solutions obtained directly 
from the Navier-Stokes equations and those obtained from 
the Reynolds-averaged Navier-Stokes equations. 
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